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The size of controllable quantum systems has grown in recent times. Therefore, the spatial 
degree of freedom becomes more and more important in experimental quantum systems. However, 
the investigation of entanglement in many-body systems mainly concentrated on the number of 
entangled particles and ignored the spatial degree of freedom, so far. As a consequence, a general 
concept together with experimentally realizable criteria have been missing to describe the spatial 
distribution of entanglement. We close this gap by introducing the concept of entanglement width 
as measure of the spatial distribution of entanglement in many-body systems. We develop criteria to 
detect the width of entanglement based solely on global observables. As a result, our entanglement 
criteria can be applied easily to many-body systems since single-particle addressing is not necessary. 

PACS numbers: 3.67.Mn, 3.56.-W, 5.50.+q 


I. INTRODUCTION 

To experimentally implement quantum technologies such as information processing, simulation [5] or metrol¬ 

ogy El: characterizing and understanding multipartite entanglement [5] is important. By characterizing multipartite 
entanglement one is able to understand experimental setups, identify experimental limitations, and investigate possi¬ 
ble noise sources in a better way. The size of controllable quantum systems has grown in recent times and large arrays 
of atoms jj>] or clouds of macroscopic singlet states [7] have been produced. Also ideas of coupling several ion traps 
to build a quantum computer exist [8j. However, for large quantum systems the spatial degree of freedom becomes 
more and more important since for large systems external fields cannot be approximated by constant fields anymore. 
Furthermore, entanglement can be easily protected again constant fields but is very vulnerable to spatially varying 
fields. As a result, the spatial distribution of entanglement compared to the spatial distribution of external fields is 
decisive for the temporal evolution of spatially extended quantum systems. 

The spatial distribution of entanglement plays also an important role in the investigation of quantum phase tran¬ 
sitions Ml and makes a distinction between different given ground states of the generalized Heisenberg spin-chain 
possible mm- 

Entanglement of multipartite system can be characterized with different quantities, such as the entanglement depth 
or fc-producibility nznm, which is defined as the minimum number of entangled particles necessary to create a 
given state. In systems without any spatial ordering, entanglement depth is a powerful variable to characterize the 
entanglement properties of this system. However, in systems with spatial order such as spin chains, or in the presence 
of gradient fields, the dynamics of a system may also depend on whether entanglement exists only between neighboring 
particles or between distant particles. 

The investigations done so far concentrated on entanglement depth or fc-producibility EMI or required address¬ 
ability of single subsystems Ml. Our criteria, developed in this paper, are based solely on global observables. 
Therefore, they open the possibility to study correlation propagation and other physical characteristics of many-body 
systems without the necessity of addressing single subsystems. 

The paper is organized as follows: First, we introduce the concept of entanglement width to characterize the 
spatial distribution of entanglement. Then, we give an example of how the width of entanglement influences the time 
evolution of a quantum system before we present methods to characterize the width of entanglement with the help of 
global observables. We conclude by demonstrating how quantum phase transitions manifest themselves in the width 
of entanglement. 


II. THE CONCEPT OF ENTANGLEMENT WIDTH 

The width of entanglement w of a pure state IT) = (£) ■ \ipj) is defined as the maximal distance w of two entangled 
particles within the states \ipj) (see Fig. [l]). A completely separable state exhibits an entanglement width of w = 1. 
The entanglement width of a mixed state is defined by the minimum with w over all decomposition g = J ~b pj |, 
that is 


w(g) 


min 

decompositions 


max {w(ipj)} 


(1) 
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FIG. 1: Comparison of entanglement depth (a) and entanglement width (b) of the state |'F) = |V'i.3,6) ® ® |^>2)- Whereas 

the entanglement depth disregards any spatial ordering, the definition of entanglement width requires the particles to be 
spatially ordered, e.g. in a spin chain. The entanglement depth of the state |'I'} in (a) is given by k = 3 (since maximally three 
particles are entangled). This is a lower bound on the entanglement width in (b), which equals w = 6 (since entanglement 
occurs over a distance of six particles in the chain). 



FIG. 2: Variance of the observable J for different parameter A and entanglement configurations for N = 16 particles. Pairs 
of encircled particles form together the state \ip~)- The product state (non-encircled particles, w = 1) is chosen in such a way, 
that it minimizes the variance. Although, both entangled states possess the same entanglement depth k = 2, they exhibit quite 
different behavior due to their different spatial distribution of entanglement. 


By definition, the entanglement depth is a lower bound of the entanglement width. However, the entanglement width 
does not make any statement about the entanglement depth. For example, the width of entanglement w = 6 in Fig. [l] 
stays the same, no matter if all particles are entangled with each other or only the two outer ones (1 & 6) whereas the 
entanglement depth changes from k = 6 to k = 2. Furthermore, states with equal entanglement depth but different 
entanglement width can lead to dramatic different effects as we demonstrate in the following example. 

Motivated by Ref. m consider a chain of particles j at the positions Xj = Xq + j ■ d, with equal spacing d between 
two particles. Furthermore, consider the observable 

N 

J = ^2 sin dj (2) 

3 =1 


with aj = 


( x) ( y) (z 

T\ ’ (J- O’- 

3 5 j i ^ J 




denoting the Pauli matrices acting on particle j and A being a parameter determining 
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the observable. Such an observable can be created e.g. by a time evolution 

t/(^)=ex (3) 

i 

of spins in a gradient field which rotates the spins depending on their position. Here, A is given by the gradient of the 
field. Another approach to create an observable similar to J is given by a standing light wave coupled to cold atoms 
in a lattice. In this way, Eckert et al. were able to distinguish between 4 given ground states with different spatial 
entanglement configuration which are important for condensed matter and high energy physics m- Whereas Eckert 
et al. needed the previous knowledge that their state is given by exactly one out of four states, the following method 
to investigate the width of entanglement does not need this strict restriction. 

J is given by the global spin operator, if d/X is an integer and e.g. xq = A/4. Its variance (AJ) 2 = (A J x ) 2 + 
(A J y ) 2 + (AJ Z ) 2 , with (AJa,) 2 = (J 2 ) — {J x ) 2 , is minimized by the global singlet state 

i*"> = 01 r)j, k (4) 

O '»*0 

were particles are combined to pairs forming together the singlet state \ip~) = (|01) — |10))/\/2. However, for d/X N 
the variance depends on the spatial distribution of the entangled pairs (j , k) as can be seen in Fig. [2] for N = 16 and 
xo = —d/2. For example for d/X = 1/(2 N) the “hugging” configuration with w = N, where particle j is entangled 
with particle N + 1 — j (solid green line), reaches the minimal variance of (AJ) 2 ug = 0. However, if all particle with 
odd number j are entangled with their right neighbor k = j + 1 (dashed blue line), the width of entanglement is given 
by w = 2, which leads to the variance 

_, Q 7r 

(AJ) r 2 n =-A(l-cos(-))«—. (5) 

As a consequence, the variance of J depends crucially on the spatial distribution of entanglement. 

Also the quantum Fisher information (QFI), which is an important measure for quantum metrology and entangle¬ 
ment [20 , is strongly influenced by the spatial distribution of entanglement. For example, the time evolution given 
in Eq. ([3 ) leads to a QFI given by F hug = N 4 for the hugging configuration, whereas it is given by F rn = 4 N 2 for the 
right neighbor configuration. As a consequence, the two here considered states exhibit different behavior although 
their are equal in their entanglement depth. 


III. FIRST CRITERION FOR ENTANGLEMENT WIDTH 


In this section we demonstrate that the observable A J, defined in Eq. ([ 2 ]), is able to distinguish between short- 
range and long-range entanglement. We will estimate the minimal variance (AJ) 2 for states with nearest-neighbor 
entanglement and demonstrate that states possessing long range entanglement are able to violate these bounds. Since 
(A J) 2 is a concave function, it reaches its minimum for pure states. Furthermore, since nearest-neighbor entanglement 
implies the entanglement of maximal two particles, we find for pure states (A J) 2 = k) w 'th the two- 

particle variance 


(AJ)fj ' k ) — 3 (a 2 + al) 


i a j(dj) + a k (d k )Y - 2aja k (djd k ) 


where aj = sm(2nXj/X). The minimum of this two-particle variance is given by 


min(A J) 2 ^ fe) = 



4e 2 

(l-s) 


S ) 2 


2 ) 


— 1 < £ < £q 
£ 0 < £ < 1 


( 6 ) 

(7) 


where we assumed with out loss of generality a.j\ > |afc| and defined e = a k /aj and = 2 — \/3 ~ 0.27 (for a proof see 
Appendix A). The remaining task is to optimize over all possible combination of entangled pairs. This is a classical 
optimization task and for many parameter A easy to estimate. For example for A = 1 /(2iV) and N = 4fc + 2 with 
t:£N the optimal pairing is given by (1, 2).., ( N — 1, N) and we find e > 2 — \/3 for all pairs of aj, a k - Therefore, the 
lower bound is exactly given by Eq. ([ 5 ]). A simple lower bound for the optimal pairing is given by 

N 1 

min (AJ) 2 > -min(AJ) 2 fe . (8) 

3 =1 
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FIG. 3: Variance of the observable J for different parameter A and entanglement configurations and N = 16 particles. Lines in 
(green, solid), (blue, dashed) and (black, dotted) correspond to the states given in Fig. [2] The limits for certain entanglement 
configurations are given by (i) (red, dashed-dotted) dower limit for the entanglement configuration (1, 2),... (15,16), (ii) (blue, 
+): lower bound for w = 2, (iii) (orange, * ): lower bound for w = 4. The state in the hugging configuration with width w = 16 
(green, solid) violates in the encircled area the lower bounds for entanglement width w = 2 and w = 4. Therefore its long range 
entanglement is detected by the presented method. 


In this way, we can find easily lower bounds not only for nearest-neighbor entanglement but also for general maximal 
widths w. For example, in Fig. [ 3 ] we show for N = 16 the variance (Aj)j( ug which can beat the limits for w = 2 and 
w = 4 for X/d « 7/(2 N). 

IV. SECOND CRITERION FOR ENTANGLEMENT WIDTH 

Another method to distinguish between short range and long range entanglement is to use several observables with 
different correlations. For example the Hamiltonian 


N 

Hi = J2djd j+1 . (9) 

3 =1 


of a spin chain in the Heisenberg model contains only nearest neighbor correlations. On the other hand, the total 
collective angular moment J c Eq. ([ 2 ]) with d/X € N includes correlations between all spins with equal weight. Both 
observables can be used to detect multipartite entanglement mmm- However, they can only together distinguish 
between nearest-neighbor and long range entanglement. 

The difference between the energy and the total spin may be a good indicator for entanglement beyond nearest- 
neighbors, which implies an entanglement width of w > 3. Indeed, by defining the correlation function 

X (N) = (J 2 c -2H)-3N 

N -2 

= EE SjSj+k (10) 

3 fc=2 


one can define a correlation measurement, which includes only non-nearest neighbor correlations. 

For quantum states with only nearest-neighbor entanglement the correlations appearing in x(N) can all be treated 
classically. The minimum of x(A) can then be derived with the help of circulant matrices |J2] (see Appendix B). As 
a consequence, all states with entanglement width w < 2 satisfy the inequality 


X (N) > -N 


sm(3n/N) 
sin(7 t/N) 


( 11 ) 
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For large N we find 


Ar sin(37r/iV) 47r 2 

— iV— 7 —-— . . ->■ —31V + . 

sm( 7 r/./V) iv»i TV 

A quantum state, which has obviously entanglement beyond nearest-neighbor entanglement is the state 

iV/4-1 

l*> nnn — (££) \lf> )fc+l,fc+31'0 )fe+2,fc+4, 

/c=0 


(12) 


(13) 


where always two non-nearest-neighbor spins form together the Bell state \4>~) = ( 101) — 110))/a/2- As a consequence, 
we find Xnnn{N) = —3 N which violates for all finite N the bound —!Vsin(37r/iV)/ sin(7r/lV). The violation decreases 
with 1/N. However, note that a maximal violation proportional to 1/N appears also in other entanglement criteria, 
see e.g. Ref. [22.. The reason may be, that our criterion is highly symmetric and the minimal overlap between 
maximal entangled symmetric states and product states decreases faster than 1/N [23], which is a consequence of the 
de Finetti theorem. 

As an application for this criterion, we consider spin chains with nearest-neighbor and next-nearest-neighbor cou¬ 
pling described by the Hamiltonian 


N N 

H 2 = J2 1 +a Yl ^ 21+2 • ( 14 ) 

l=i l=i 


This model is often called the Ji-J 2 -model. It is used to understand phenomena in magnetic materials such as 
Tomonaga-Luttinger-liquids states and spin-Peierls states [251 :2Sj. Here, the interactions between nearest and next- 
nearest neighbor compete with each other, leading to frustrated spins. 

To estimate the quality of our criterion, we compare it with a criterion for entanglement depth. For separable state 
the minimal energy ( H 2 ) sep can be computed with the help of the eigenvalues of the circulant correlation matrix. 
With this method, we determine for separable states (H 2 ) sep > -^fcdrcuiant with 


hcircnlant — mill 
m 


TYl /Til 

cos( 27 r—) + acos( 27 r-^-) 


(15) 


To estimate the minimal achievable energy H 2 for states with entanglement depth fc < 2 we use the methods from 
Refs. IB’, m ( see Appendix C). As a consequence, for states with entanglement depth k < 2 we find the limit 
{H 2 ) 2 -p ro d / Nh 2 p ro( \ with 


^2prod — 1 T Cl H~ 


1 

2 +4a' 


(16) 


As an example we investigate the entanglement of the ground state of a spin chain with N = 8 spins and an 
interaction Hamilton given by H 2 defined in Eq. (14) for different values of a [32]. 

In Fig. [4] we compare the expectation value of the energy H 2 (blue crosses) with the minimal achievable energy 
for states with entanglement depth k < 1 (blue dashed dotted line) and k < 2 (blue dashed line) as well as the 
correlation function x (red *) and the limit for states with entanglement width w < 2 (red line). As can be seen in 
Fig. 0 the ground state exhibits entanglement depth of fc > 3 for a < 0.3 and a > 1.1. On the other hand, we detect 
entanglement width of w > 3 for 0.86 < a using y. As a consequence, there exist a regime (black square) where 
we detect already an entanglement width of w > 2 with the help of the correlation function \i but no multipartite 
entanglement with existing methods. Furthermore, the correlation function x shows two sudden jumps at a ~ 0.5 
and a « 0.7 which are indicators that for N —> 00 phase transitions may occur at these points [12] . Contrary to that, 
we find no hint for phase transitions in the expectation value (H 2 ). 

The behavior found with our entanglement criteria fits well to existing results in the literature. As shown by 
Majumdar and Ghosh [27[ [28], the ground state of H > is given by 


V!>mg) 


\i> ) 


2k,2k±l 


(17) 


for a = 1/2. Therefore, the ground state at this point shows bipartite entanglement between nearest neighbors as also 
indicated by our results. Furthermore, the relation between entanglement and phase transitions for the Hamiltonian 
H 2 has been investigated with the help of the concurrence [12] and generalized geometric measures m- Similar to 
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FIG. 4: Comparison of entanglement depth and entanglement width for different a and N = 8 spins: the energy H 2 of the 
ground state (blue “+”) indicates entanglement if it lies below /i c i rC uiant (blue dashed-dotted line) and multipartite entanglement 
if it lies below foprod (blue dashed line). The expectation value of \2 (red “*”) indicates entanglement beyond nearest neighbor 
if it lies below Xciass (red line). As a consequence, our results detect entanglement of non-neighboring particle before we detect 
multipartite entanglement (black square). 


our results, Gu et al. found with the help of the concurrence that the entanglement changes from nearest neighbor 
to next-nearest neighbor entanglement around a ~ 0.75 for N = 8 spins jT2]. Biswas et al. showed that phase 
transitions which occur only for N —> 00 manifest themselves in the generalized geometric measure for entanglement 
also for finite N around a ss 0.7 m- Our entanglement criterion is easily accessible experimentally in contrast to 
the methods used by Biswas et al. and Gu et al., which need the exact knowledge of the state and addressability of 
single particles. 


V. CONCLUSION 

In summary, we have introduced the concept of entanglement width and demonstrated that states with equal 
entanglement depth but different entanglement width behave differently. We developed criteria based solely on global 
observables which are able to distinguish between different values of entanglement width. With the help of these 
entanglement criteria we investigated the ground state of a spin chain described by H 2 and showed that we are able 
to detect long range entanglement before we detected multipartite entanglement. 

Furthermore, phase transitions as predicted by other theories for spin-chains described by H 2 manifest themselves 
in our concept of entanglement width whereas they were not visible in our measure for entanglement depth. Although, 
we just started to investigate the width of entanglement it turned out to be an important indicator to investigate 
the behavior of quantum systems in the presence of space dependent interactions. In this way, we could used it as 
indicator for many-body phenomena such as phase transitions. It would be very interesting to further investigate this 
possibility. 

In addition, detecting the width of entanglement with global observables is not limited to the two methods demon¬ 
strated in this paper. Another example is e.g. given by measuring the variance of the total spin with and without a 
linear dependency on the position and a consecutive comparison. Here, the variance of the total spin without exter¬ 
nal field gives information about the entanglement, whereas the position dependent observable measures the spatial 
distribution of correlations, classical as well as quantum (for more details see Appendix D). We are certain that many 
more criteria detecting the width of entanglement will be developed in the future. 

We thank G. Toth, R. Sewell, J. Kong and M. Mitchell for fruitful discussions. This work has been supported by 
the EU (Marie Curie CIG 293993/ENFOQI), the FQXi Fund (Silicon Valley Community Foundation), the DFG and 
the ERC (Consolidator Grant 683107/TempoQ). 
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Appendix A: The minimal two-particle variance 


The two particle variance (A J) 2 . k ^ = 1 + e 2 — g(e, \ip)) with 

9 (e, IV’)) = ((ctj) + £(CTfe)) 2 - 2e(a j a k ) (Al) 

is minimal if g riches its maximum. To maximize g we parametrizes the state \ijj) by 

| ip) = alV’ - ) + 6|i/> + ) + c\4>~) + d\(j) + ). (A2) 

with the usual definition of the Bell states l^ - ), \'4 >+ ), \4 >~) and |</> + ). As a consequence, the spin vector of the single 
particle is give by 

/ Re(M*) TRe(ac*) \ 

(a)j/ k = 2 I Im (6*c) ± Im (ad,*) I . (A3) 

y Re (c*d) ± Re (a6*) J 


Since the two vectors always lie in a plane, we assume (a v ) = 0 for both particles. This implies phase relations 
between the parameters a — d which implies that the maximum of g can be reached by choosing all parameters to be 
real. As a consequence we get 

g(e, \ip)) = 4(6 2 + c 2 )d 2 ( 1 + e) 2 + 4(6 2 + c 2 )a 2 (l — e) 2 + 8 ea 2 — 2e (A4) 


which does not depend on b 2 or c 2 itself but only on their sum. Therefore, we choose w.o.l.g 6 = 0 and arrive at 


g(e, | ip)) = 8a 2 e(l — 2c 2 ) — 4c 4 (l + e) 2 — 2e 


(AS) 


where we also made use of the normalization relation a 2 + b 2 +c 2 + d 2 = 1 of the state. The function g is positive and 
monotonically increasing in a 2 for e > 0 and c 2 < 1/2 and negative for £ > 0 and c 2 > 1/2. Therefore, the maximum 
of g for e > 0 is reached by choosing the maximal possible value of a 2 given by a 2 = 1 — c 2 . Now, we are able to 
maximize g over the only left parameter c 2 . The maximum is reached for 


c 


2 


1 £ 

2 “ (I^F 


(A6) 


However, this is only possible for e < 2 — a/ 3 since c 2 > 0 . If £ > 2 — v/3, we have to choose c 2 = 0 and a 2 = 1 to 
maximize g. As a consequence, for e > 2 — \/3 « 0.27 the variance A J is minimized by the singlet state \tp~)- 

In a si milar way, the maximum of g for £ < 0 can be estimated and we find that choosing again c as defined in 
Eq. (A6) is optimal. However, for —1 < £ < 0 we find 1/2 < c 2 < 3/4 which makes now further distinction of cases 
necessary. 


Appendix B: Lower bound of y(lV) for w = 2 


For w = 2 the correlation function x(7V) defined in Eq.(5) becomes classical. Therefore, it can be written as 
X cl {N) = X T CX, with X = (xi,... ,xn) t , Xj = (a) and the circulant correlation matrix C. As a result, the 
minimum of Xjv f° r only nearest-neighbor entangled states is bounded from below by the minimal eigenvalue of C. 
The correlation matrix C is a circulant matrix [16], and the eigenvalues of C are given by 


A 


m 



n 


N - 3 

sin(TT^) 


for m = 0 

for to = 1,..., N — 1 


(Bl) 


with c n denoting the entries of the circulant matrix C. As a consequence, the minimal eigenvalue of C is given by 
min(A m ) = — sin(37r/iV)/sin(7r/lV). 

To estimate the minimum of X cl (A0 we have to consider also the length of the eigenvector X m . Since Xj < 1 the 
maximal length of the vector X is given by N. Indeed, by choosing X = (a?i,... ,Xn) t with e.g. 


/ cos2t\ 
sin 2-7T 

V o J 


(B2) 





we have found an eigenvector of C with the minimal eigenvalue and the length N. 

We note two interesting facts: (i) the eigenvectors of circulant matrices like C are independent of the coefficient 
Cfc, only the eigenvalues depend on c*,. (ii) Although the vectors Xj live in a three-dimensional space, the vectors Xj 
forming the eigenvector are only two-dimensional. Since the eigenvectors are independent of the coefficient c& the 
scheme presented here can be generalized to all correlation functions 

x = £E Cj,kXjXk (R3) 

3 k 

which are circulant, that is Cj^ = cu_u and the minimum of x is always given by N ■ min(A). As a consequence, the 
scheme presented here can be used to find the minimal energy of arbitrary circulant Hamiltonians. 


Appendix C: Minimal energy H 2 for k = 2 


To derive the minimal Energy (H 2 ) achievable by states with entanglement depth k < 2 we divide the spin chain 
into blocks, where the ions within a block may be entangled but ions belonging to different blocks must be separable 
analogue to Refs. [11,12]. Using the methods of Ref. [12], the minimal achievable energy is bounded by 

(ci) 

3 

where Cj describes the optimization over each single block. For single-particle blocks Cj is given by 

^single _ max. (x 2 + axl ) = 1 + a. (C2) 

For two-particle blocks we obtain with the help of Lemma Al of Ref. [12] 

-x k x t + -(a?fc +xf) + a(xl + xf) 

= 2(1 + a) + TT2W (C3) 

As a consequence, for states with entanglement depth fc < 2 we derive the limit (l? 2 ) 2 -prod > —Nh 2 plo d with 




double 


= max 

lb) 


/l 2 prod - 1 + a + 2 + 4c U 


(C4) 


Appendix D: Detecting the width of entanglement with the help of a gradient 

Another way to detect the entanglement width is to change the weight of different interactions through an external 
space-dependent field, e.g. by measuring the collective spin in the presence of a magnetic gradient. In the following, 
we assume for simplicity an even number of spins and that the magnetic field is zero in the middle of the chain. 
Therefore, the interaction of the magnetic field with particle k can be described by 

BjP = {2k-N-l)ajP , (Dl) 

with £ = x,y, z. Similar to the angular momentum we define and B 2 = (B x ) 2 +(B y ) 2 +(B z ) 2 . Whereas 

small values of (J 2 ) are indicators for entanglement, large values of (B 2 ) are indicators for long range correlations, 
both classical and quantum. As a consequence, all four-qubit states which are separable under the partition 1,213,4 
obey the inequality 


(B 2 ) < 24 + 10(J 2 ). (D2) 

as depicted in Fig. [ 5 ] This bound is tight, since it is reached at the point (( J 2 ) = 0, (B 2 ) = 24) by the state 
|^ _ )i,2|V’ _ )3,4 and at the point ((J 2 ) = 8 , (B 2 ) = 104) by the state | tt)i,2| 44 -)3,4- Mixtures of these two states lie 
exactly on the line. 
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FIG. 5: Expectation values of J 2 and B 2 for states of the form: |'i/’}i,2|V’)3,4 (•), |V')i,2|</ > }3,4 (•), and |t^)i,4|^)2,3 (□). States 
which are separable under the partition 1,213,4 lie below the limit: {B 2 ) =24+ 10(J 2 } (red line) 



Eq. (D2| can be proven with the help of semi-definite programming. Here, we utilize the fact that all states which 


are separable under a gi ven partition possess a positive partial transposition (PPT) under this partition [23 EB- 
Therefore, we prove Eq. (D2) by searching the minimum of 10(J 2 ) — ( B 2 ) for four-qubit states g under the condition 
that the partial transpose of g under the partition 1,2|3,4 is positive semidefinite. 

However, also states which are separable under the partition 1|2,3|4 possess only nearest-neighbor entanglement, 
too. For these states, an inequality similar to Eq. (D2) can be formulated. Here, we search the minimum —a of 
m(J 2 ) — (B 2 ) for four-qubit states g which are PPT under the partitions 112, 3,4 as well as 1,2,3|4 for different 
to. In this way we get a family of inequalities (B 2 ) < a + m ■ (J 2 ). In contrast to Eq. (D2), there exist no optimal 


inequality: either we optimize the y-intercept a or the slope to. Furthermore, we made a strong relaxation by going 
from separability under the partition 112, 3|4 to PPT under the partitions 112, 3,4 and 1, 2, 3|4. Therefore, the gained 
inequalities are not tight. Nevertheless, we find joined inequalities for states which are PPT under the partition 
1, 2|3,4 or simultaneous PPT under 112,3,4 a nd 1 , 2, 3|4 for example (B 2 ) <24+ 16.1 • (J 2 ) (optimal y-intercept) or 
( B 2 ) < 50.2 + 10 ■ (J 2 ) (same slope as in Eq. (D2|). All states which violate at least one of these inequalities possess 
entanglement beyond nearest neighbors. 

Similar bounds can be found for TV = 6. For example we find the bound (B 2 ) < 36 + 48(J 2 ) for states which 
are simultaneous PPT under the partition 1, 2|3,4, 5, 6 and 1,2, 3,4|5,6. Unfortunately, such optimization problems 
are very complicated for large systems and can therefore not be solve with semi-definite programming for large TV. 
However, for (J 2 ) = 0 and (J 2 ) = 2TV we are able to analytically calculate the upper bound of (B 2 ), leading to a 
conjecture about the general behavior: 

The only states with (J 2 ) = 0 and only next-neighbor entanglement are given by |T _ ) = (££) \'4’~)2j-i.2j which 

leads to (H 2 )_ = 6TV. The state |^ cl ) = | f, ■ ■ •,t>i .JV/ 2 I I, • • ■ ,1)n/2+i,.. . )JV possess the highest possible value of 

(B 2 ) for states with only entanglement within the subgroups {1,..., iV/2} and {TV/2 + 1,...,TV}. This state is 
characterized by the expectation values (J 2 ) c i = 2TV and (B 2 ) c i = (TV 4 )/4 + 2(TV — 1)TV(TV + l)/3. Therefore, all 
states with expectation values ( B 2 ) exceeding this bound must posses entanglement between these two subgroups. 
If we allow entanglement only between nearest neighbors, then compared to (B 2 ) c 1 only the correlation between 
the two middle particles within the chain are allowed to change from two to six. As a consequence, all states with 
(B 2 ) > (B 2 ) c i + 4 posses entanglement beyond nearest neighbors. 

Furthermore, the behavior of (B 2 ) for TV = 4 leads us to the conjecture, that the largest possible value of (B 2 ) for 
a given (J 2 ) with only entanglement within within the subgroups {1,..., TV/2} and {TV/2 + 1,..., TV} is reach by the 
state g = p|4' _ )(4'”| + (1 — p)|4' cl )(4' cl |. As a result, we obtain the conjecture that all states with entanglement only 
within the subgroups {1,..., TV/2} and {TV/2 + 1,..., TV} satisfy the inequality 


(B 2 ) < (H 2 )_ + 


(B) cl - (B 2 )_ 
2TV 


(A 


(D3) 


For TV = 4 the conjecture coincides with the inequality given in Eq. (D2). Furthermore, first tests with random states 
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confirm this bound. 
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